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Although a great deal is known about design techniques for optimum 
(in a minimax error sense) finite impulse response (FIR) low-pass digital 
filters, there have not been established any practical design rules for such 
filters. Thus, a user is unable to easily decide on the (approximate or exact) 
filter order required to meet his design specifications and must resort to 
tables or trial and error procedures. In this paper, such a set of design 
rules is given. In the case of very narrow bandwidth or very wide bandwidth 
filters, analytic relations between the filter parameters can be readily ob- 
tained. In all other cases, exceedingly good linear and nonlinear fits to the 
data can be obtained over somewhat restricted ranges of the parameters. 
These fitting procedures lead to a practical set of simple design rules for 
estimating filter order from the desired specifications. 

I. INTRODUCTION 

The problem of designing an optimal (in the minimax sense) low- 
pass FIR digital filter to meet design specifications has been thoroughly 
investigated 1-5 and may be considered to be solved. Thus, given a 
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Fig. 1 — Definition of low-pass filter parameters. 

specified impulse response duration of N samples, a passband cutoff 
frequency F p (see Fig. 1), a stopband cutoff frequency F„ and a ratio 
of tolerances K = Si/ 82 (where 81 = passband tolerance and 82 = stop- 
band tolerance), an optimal approximation to these specifications can 
be designed. The approximation is optimal in the sense that, for given 
values of N, F p , F„ and K, 81 (or equivalently 82) is minimum. The 
nature of the solution is such that there are three distinct classes into 
which it may belong, depending on the specific design parameters. 
These classes have been called extraripple filters, 4 ' 6 scaled extraripple 
filters, 7 and equiripple filters with one less than the maximum possible 
number of ripples. The differences between these classes lie in the 
number and amplitude of the ripples in the weighted error curve. The 
weighted error curve is defined as: 



1 - H(e**f) 



E(e* Tf ) = 
where H(e j2 ' f ) is the frequency response of the optimal filter. Extra- 



K 

-H(e* rf ) 



^ / ^ F p 
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ripple filters have (N + 5)/2 equal magnitude extrema in their error 
curves. Scaled extraripple filters also have (iV -j- 5)/2 extrema, all but 
one of which are of equal magnitude. The third class of solution has 
(N + 3)/2 equal amplitude extrema in its error curves. Figure 2 shows 
plots of curves of transition bandwidth, AF = F s — F p , versus pass- 
band cutoff frequency, F p , for two sets of conditions. 6 The data in Fig. 
2a show the curves for K = 1, S t = S 2 = 0.1, JV* ■ 9, and AT =11; 
whereas Fig. 2b shows the curves for K = 100, 8 t = 0.01 (8 2 = 0.0001), 
N = 19, and AT = 21. The minima along each curve are the extraripple 
filters. (For fixed values of <5i and 8 2 , as in this figure, there are only 
(N — l)/2 distinct extraripple filters.) The local regions around the 
minima are the scaled extraripple filters 7 and the remainder of the curve 
represents equiripple solutions with one less than the maximum number 
of ripples. As will be shown in the next section, the first and last extra- 
ripple filters can be obtained analytically because they are simply re- 
lated to the Chebyshev polynomial of appropriate degree. 89 

As seen from Fig. 2, as F p varies from to its maximum possible 
value, the transition width goes through a sequence of minima and 
maxima. The variation in the transition widths of the minima and 
maxima decreases as N increases. Furthermore, the variation in transi- 
tion width between adjacent maxima and minima also decreases as N 
increases. In fact, except for a narrow region at the beginning and end 
of the curve, the curve of transition width versus passband cutoff fre- 
quency is relatively flat over a wide range of values of F p , Si, and K. 
Figure 3 shows a sequence of three plots of transition width versus 
passband cutoff frequencj^ for extraripple filters (i.e., only the minima of 
the curve are plotted) of length N = 101 for various values of 8i and 
K. (The entire curves are not plotted because the amount of computa- 
tion required for a smooth curve of such high order is inordinately 
high.) Figure 3a shows a sequence of four curves for K = 1, 5i = 0.1, 
0.01, 0.001, 0.0001, whereas Fig. 3b shows the same sequence for 
K = 10, and Fig. 3c shows the sequence for K = 100. Several observa- 
tions can be made from this figure. 

(i) For a wide range of values of F p , and fixed 8i and 8 2 , the transi- 
tion width of the extraripple filters is relatively insensitive to 
F p . The larger the value of K, or the smaller the value of 5i, 
the worse this approximation becomes. 

(it) In the regions of either very small or very large values of F p , 
the transition width generally decreases. 

Based on these curves and on previous results with window designs, 10 
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Fig. 2 — The transition width as a function of passband cutoff frequency for two 
sets of filter parameters. 
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Fig. 3 — The transition width as a function of passband cutoff frequency for 
N = 101-point extraripple filters and K = 1, 10, and 100. 

it seems fairly reasonable to expect some simple design relationships 
to exist between the five basic filter parameters, N, F p , F B , 8 h and <5 2 
(or K), at least in the extreme case of the Chebyshev solution, and also 
for a reasonably large region near F p = 0.25. Experience in practical 
situations has shown that the number of terms needed in the optimum 
FIR low-pass filter to meet design specifications is significantly less 
than the number of terms estimated by known relationships on win- 
dows. 10 Therefore, the goals of this paper are to obtain approximations 
to the actual design relationships between linear-phase, low-pass filter 
design parameters and to illustrate their use in actual design examples. 
The organization of this paper is as follows. First, the design rela- 
tionships for the Chebyshev solution are derived, and approximate 
formulas for the transition width in the limit of large values of N are 
obtained. Then the results of measurements on a wide range of filters 



774 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1973 

where the number of passband and stopband ripples are equal are 
shown. Then minimum mean-square relative error fits to the observed 
data (for large N) assuming both linear and nonlinear dependency on 
the basic filter parameters are computed. To apply the design relation- 
ships for all values of N, a correction formula is derived, based on con- 
sideration of the appropriate transition width of the filter. Finally, 
a set of rules is presented for going from a set of desired filter param- 
eters to an estimate of the lowest-order filter which meets these 
specifications. 

1.1 Summary of Design Relationships 

Given the low-pass filter parameters F p , F„, 81, and 8 2 (or, equiva- 
lent^, K = 81/82), the minimum filter impulse response duration, N 
required to meet the above specifications can be estimated from the 
relation 



„ DJM, fr) - f(K)(F s - F,)* 

" " " (F B - F~) ~ ~ + L ' 



where 



D m (8 u 8 2 ) = [5.309 X 10- 3 (logio 5i) 2 + 7.114 X 10" 2 (logio 5,) 
- 0.4761] logio 8 2 - [2.66 X 10" 3 (logio 5i) 2 

+ 0.5941 (logio 81) + 0.4278] 

and 

f(K) = 0.51244 logio K + 11.01217. 

The above relations are valid to within 1.3 percent relative error in N 
if 81 £ 0.1 and 82 ^ 0.1. 

II. CHEBYSHEV SOLUTIONS 

Let \h(n), n = -(N - l)/2, • • • , (N - l)/2} be the impulse re- 
sponse of the desired digital filter. (N is assumed to be odd throughout 
this paper.) The impulse response satisfies the symmetry condition 
h(n) = h(—n) to give the desired linear phase. The frequency response 
of the filter is given by 

H(e**') = M0) + E 2h(n) cos (2t fn). (2) 

n=l 

By making the substitution 

cos (2ir/) = ,^ , h ( 3a ) 



m 



OPTIMUM FIR LOW-PASS DIGITAL FILTERS 775 

or 

^ = (^ ti )cos(2 I r/) + (^ =J ) (3b) 

the interval ^ / ^ 0.5 is mapped to the interval X ^ x ^ —Lit 
is easily shown 37 that the mapping transforms the trigonometric poly- 
nomial of eq. (1) to an algebraic polynomial in x of the form 

(-V-D/2 

G(x) = £ b(n)x», (4) 



i=Q 



where the \b(n)\ are straightforwardly related to the \h(?i)\. 

The basic filter design problem is to find coefficients b(n) [or fc(n)] 
such that the weighted error of approximation is equiripple in both the 
passband and stopband. In the case when either the passband or the 
stopband has only one ripple, the solution to the filter design problem 
may be obtained analytically based on the theory of Chebyshev poly- 
nomials. In all other cases, alternative techniques must be used to ar- 
rive at the appropriate solution. 

Consider the Mth. degree Chebyshev polynomial, T M (x), as shown in 
Fig. 4 for M = 4. The standard representation for T M (x) is 

T M (x) = cos \_M cos _1 :r], \x\ ^ 1 (5) 

and 

T M (x) = cosh [ilf cosh" 1 sj, \x\ > 1. (6) 

In the interval — 1 £ a; £ 1, T M (x) oscillates between ±1 and, be- 
yond the value x = + 1, T M (x) grows approximately as x M . If we define 
X (see Fig. 4) as the point where T M (X ) = (1 + 5i)/8 2 , and X p as 
the point where T M (X P ) = (1 — 5i)/8 2 , it is readily seen that 8 2 - T M (x) 
is a polynomial of the form of eq. (4) [with M = (JV — l)/2] which 
is equiripple in both the passband X p ^ x ^ +X , and the stopband 
-1 S i S +1| and hence is an optimal solution to the filter design 
problem. Of course, this solution is a special case of the general solution 
in that there is only one ripple in the passband, but at least it is an 
analytically tractable case from which a great deal can be learned about 
the relationships between the filter parameters. 

It is straightforward to solve for the points X p and X in terms of 
Si, 8 2 , and N. If we set M = (N - l)/2, then at x = X > 1.0, we get 
the relation 

T M (X ) = ^4^ - cosh [M cosh" 1 Z ] (7) 

02 

or 

Z = cosh T jj cosh" 1 ( 1 -±1± \ 1 . (8 ) 
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Fig. 4 — A fourth-degree Chebyshev polynomial denning the points X<, and X v . 

At x = X p > 1.0, we get the relation 

1 - 8 



or 



T M (X P ) = ^-r-2 = cosh [AT cosh- 1 X p ~] 

02 

X p = cosh [ ^ cosh- (^-^)j- 



(9) 



(10) 



The inverse mapping of eq. (3) can be used to determine the filter 
cutoff frequencies by the relation 



, 1 _, \2x -X,+ 1 1 
/= 2^ C0S [ Zo+1 J' 



\x\ ^ X . 



(11) 



Thus, F p (corresponding to X p ) and F, (corresponding to x = +1.0) 
can be readily obtained using eq. (11) with the appropriate values for 
x. In this case, the transition width (F. - F p ) can be analytically de- 
termined for all values of N, 8 h and 8 2 as: 



"-'•-'•-sMJrri) 



— cos 



-'( "'r/i. +1 )]- (i2) 
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In the limit of large values of N (or, equivalently, M), eq. (12) can be 
simplified by the following approximations. First, for sufficiently 
small values of a, 

a- 



cosh a ~ 1 + 



2 ' 



a« 1. 



(13) 



Thus, eqs. (10) and (8) can be simplified to the form [replacing M by 
(N - l)/2] 



x, 



X, 



1 + 



jo^iH-mr 



The approximation is then made that for sufficiently small e, 

cos -1 (1 - e) « V2l. 
Thus, F s and F p are well approximated as 

cosh" 



2tt \ N - 1 



'.-£(^)[[-~mi 



(14) 



(15) 



(16) 



(17) 



cosh 



i^MT (18) 



Thus, for large A'' (N ^> 1) the approximation to the transition width 
curve is given by 

1 



AF = F s - F, 



[[ 



— cosh 



ir(N - 1) 

1 + «i 



)r-K i ( L ^)]T- ™ 



(This approximation is valid to within 1 percent for N ^ 51.) Equation 
(19) shows AF to be inversely proportional to (N — 1). This identical 
inverse behavior has been noted previously for filters designed by 
windowing. 10 These and other considerations lead one to consider as 
a performance measure of a low-pass digital filter the quantity D de- 
fined as 

D = (N - 1)AF = (N - 1)(F, - F p ) (20) 

which, in many cases, depends only on 61 and 5 2 . 
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The curves in Fig. 5 (for JV = 127) show plots of performance D 
versus 20 Iogio (8 2 ) both for values of K = 5i/5 2 from K = 1 to 
K = 1000 and for values of 8i from 0.5 to 0.0001. The behavior of these 
curves is as predicted from eq. (19) and from intuition about the be- 
havior of D as 5i and 8 2 get large. It is clear that when 5i + <5 2 = 1.0 
[i.e., 8i = K/(K + 1), 8 2 = 1/(K + 1)], then D = since there is no 
transition band. In this case, the term cosh -1 (1 — 8i)/8 2 vanishes, and 
the first and second terms in eq. (19) cancel exactly. In the limit of 
small values of Si, the second and third terms in eq. (19) are approxi- 
mately equal and cancel. Thus, D is approximately of the form 

D«icoBh-»(l+! x )- (21) 

Since 

cosh" 1 (x) = In (x + Vz 2 - 1), (22) 

eq. (21) can be rewritten as (assuming Si negligible) 

D « - (In 2 - In 8 2 ) (23) 

which is independent of St. Thus, as seen in Fig. 5b, in the case of small 
Si, D is essentially independent of Si. 

The curves of Fig. 6 show the behavior of D for various values of 
JV ^ 127, for K = 1, 10, and 100. The values of JV used for these plots 
were N = 3, 7, 11, 19, 51, and 127. The differences between the data 
for N = 127 and the data for N = 51 are relatively small. These 
curves also exhibit another interesting phenomenon. The curve for 
N = 3 saturates at a value of D = 1. This is due to the limitation that 
the transition width, A.F, must be less than or equal to 0.5. Thus, the 
saturation value of D is (JV — l)/2, or 1.0 for JV = 3. The larger the 
value of K, the larger the value of 8 2 beyond which the curve for JV = 3 
saturates. 

In the case of the Chebyshev solution to the optimal filter design 
problem, a formula can be derived for the impulse response duration 
JV of a filter whose response meets specified values of Si, 8 2) and F,. 
Since F p is not specified, this result is useful only as a first guess of 
a value of JV which meets specifications on all four filter parameters. 
From the discussion given earlier in this section, / = F 8 when x = 1.0. 
Thus, eq. (11) can be used to solve for X as 

_ 3 - cos (27nF 8 ) , . 

Xo ~ H-oob(2t/?.)' {ZV 
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Fig. 6 — Plots of D versus logio 5 2 as a function of N for K = 1, 10, and 100, for 
the Chebyshev solutions. 
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Equation (8) may now be used to solve for N as 

-7T3 

L 1 + cos (2a-F.) 
Using trigonometric identities, eq. (25) can be simplified to the form 

cosh -.(L+M 

N = 1 + 7i\ (26) 

cosh ~' (c"Sw ) 

For these values of X and 2V, F p may be obtained from eqs. (10) and 
(11). If the transition width obtained is too large (too small), AT is then 
decreased (increased) until the desired specifications are approximately 
achieved. (Since there is only one Chebyshev solution for fixed Si, 5 2 , 
and N, exact values of both F p and AF cannot usually be obtained.) 
As will be seen later, eq. (26) forms the basis for estimating a lower 
bound on the filter order required to meet given design specifications 
on AF, Si, and 5 2 . 

As mentioned earlier, Chebyshev polynomials can be used as the 
optimal solution in the case when there is one ripple in either the pass- 
band or the stopband. To see how this second case can be handled, con- 
sider the Chebyshev polynomial T M (x) where T M (X ) = (1 + $ 2 )/$i 
and M = (N — l)/2. As shown earlier, this polynomial represents an 
optimal filter with passband ripple S 2 , stopband ripple $i, passband 
cutoff frequency F p , and stopband cutoff frequency F,. Consider the 
function Rm{x) defined as 

R M (x) = 1 - tiT M (-x + X - 1), (27) 

where — 1 ^ x ^ Xo. An examination of the properties of Rm(x) 
shows 

(i) Rm(x) is a polynomial in x of degree M. 
(ii) In the interval X — 2 ^ x ^ X , R M (x) oscillates between the 

values 1 — $i, and 1 + $ t , 
(in) In the interval -1 ^ x ^ X - 1 - X p , R M (x) goes from 

— 52 tO 02. 

If we define F p and P. as the equivalent filter cutoff frequencies, then 
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it is readily shown that 

P, = 0.5 - F a 
F, = 0.5 - F p . 

Thus R M (x) is a polynomial with only one stopband ripple which satis- 
fies the filter optimality criterion. In summary, to design Chebyshev 
approximations with only one stopband ripple, one merely makes the 
substitution 

8i = 02 

52 = Sl „ (29) 

F p = 0.5 - F. 

F. = 0.5 - F p 

and solves for an equivalent filter with one passband ripple using the 
formulas of this section. 

The data on the Chebyshev solutions provide valuable insights into 
the behavior of the design relationships between filter parameters in 
more general cases. These data and their design relationships are pre- 
sented in the next sections. 

III. MEASUREMENTS OF D 

Earlier, it was shown that the transition width or, equivalently, the 
performance measure D for the Chebyshev solutions (i.e., either small 
F p or large F„) was generally significantly smaller than for most of the 
range of values of passband cutoff frequency. To obtain data on a more 
realistic set of filters, the value of D was measured for a large number of 
extraripple filters designed with the constraint that the number of pass- 
band and stopband ripples were the same. In this manner, the passband 
and stopband widths were almost equal and, thus, the measured data 
would characterize the parameter set over as wide a range of values of 
5i, 5 2 , and F p as possible. Six different values of N were used including 
JV = 3, 7, 11, 19, 51, and 127. Over 1500 filters were designed to 
cover the parameter range 0.00001 ^ 5 2 ^ 0.5, and 1 ^ K g 500 
(K = Si/ 8 2 ). Figures 7a through 7i show plots of D versus 20 logio (5 2 ) 
for the nine values of K and the chosen values of N. All these data are 
presented since they are fairly general and may be useful in a wide 
variety of contexts other than this paper. (Also, their measurement re- 
quired almost 3 hours of computer time on a fairly fast processor.) 

It is remarkable how similar the plots of D versus log 8 2 for the more 
general case of Fig. 7 are to the identical plots for the Chebyshev solu- 
tions. Similar behavior for small D is expected, since D tends to as 
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Si 4- 8 2 tends to 1.0 independent of N, F p , and F 8 . However, the ap- 
proximately linear behavior of D as a function of log 8 2 and log K (for 
large N) is unexpected. Another similarity between the two cases is the 
independence of D of (N — 1) for large N. The tendency of D to 
saturate for small values of N is yet another similarity between the 
curves. The main difference between the sets of curves is that, in the 
Chebyshev case, D is approximately independent of 81 for small 8y. 
This behavior is not observed for the extraripple filters of Fig. 7. 

A summary of the behavior of D for N = 127, as a function of log 8 2 
for various values of either K or Si, is presented in Fig. 8. The values of 
5i used were 0.5, 0.2, 0.1, ■ • • , 0.00002, whereas K ranged from 1 to 
500 as in Figs. 7a through 7i. In some sense, this figure represents a set 
of design curves for high-degree low-pass filters. In the next section, 
we show how the data of Fig. 8 can be approximated by linear and non- 
linear fits, and how simple modifications can be made to correct the 
results for values of N less than 127. 

IV. DATA-FITTING PROCEDURES 

In order to make most efficient use of the data of the previous sec- 
tion in a practical design problem, it is useful to express the relation- 
ships between the filter parameters in a simple manner. Since we know 
of no way of deriving exact analytical formulas, as in the Chebyshev 
case, a minimum mean-square relative error fit to the data over a re- 
stricted but reasonable range was sought. Both a linear and a nonlinear 
fit to the data of Fig. 8 (N large) were obtained. Corrections for smaller 
values of N were then obtained giving a complete set of design rules. 

The data of Fig. 8 suggest that except in the region D = (large 
values of 8i) a simple linear fit can be obtained. The curves of Fig. 8 
were assumed to be of the form 

■Olin = a -f b logio (<5 2 ) + c log™ K, 
a = -0.803, 

b =-1.359, <*» 

c = -0.737. 

where Dlin is the predicted value of D. The values for a, b, and c were 
chosen to minimize the sum of the squares of the relative differences 
between Dlin and D for the N = 127 data for values of <5i in the range 
0.01 ^ Si; 0.01 ^ 5 2 . The reason relative rather than absolute errors 
were considered is that a fixed percentage error in D, 8D, approximately 
gives a fixed percentage error in N, SN, when transition width is held 



784 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1973 




-100 -90 



-80 



-70 



-60 -50 -40 

20LOG, (fi 2 ) 



-30 



-20 



-10 



Fig, 7— plots of D versus logio S 2 as a function of N f or K = 1, 2, 5, 10, 20, 50, 100, 
200, 500 for the case of extraripple niters with the same number of passband and 
stopband ripples. 



OPTIMUM FIR LOW-PASS DIGITAL FILTERS 



785 



^ 






(e) 


- 


^"""77"- — >-^^^^_ 


N p =N 5 








^^ K = 20 
1 I^^V^ 1 


i 




3 " ■ 

1 1 1 1 





N^** 

^^L 


N p =N s 


(f) 




-— Z*^^^*s^ 


K = 50 




— 


rr~- — --^^^^^. 






- 


"?" " — ^^-^T^^^s^ 


1 ^^. 1 1 


1 




3 " ■ 

1 1 1 1 



3 - 



2 - 



o 
li- 
ar 1 







N p = N s 


(g) 


— 


iT — ^^^. 


K = 100 






7* — — -—Z^^^^^- 


T^^j i i 


" 




3 ~ — __ 

1 1 1 1 









(h) 




"""-"^T^^^ 






- 


7~" "" ^--"^^^^ 


N p = N s 
K = 200 






3 _^_^ 
1 1 1 1 


"**^^ 1 1 1 


i 







(i) 




"^^^^^s. 




^*- 


-.^Tr^C^^^ 


N p = N 5 




7"- -"-^o^^^ 


K = 500 




3 — """"*— —-»J^ 

i i i n^ 


i \ 1 1 1 1 



-100 -90 -80 



-70 -60 -50 -40 

20LOG 10 (6 2 ) 

Fig. 7 (continued). 



-30 



-20 



-10 



786 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1973 

6 




Pi gi 8— Plot of D versus logio 5 2 for the N = 127 data of Fig. 7 showing curves for 
various values of K and Si. 



fixed, as in most design problems, i.e., 



8D 
D 



8(N - 1) 
N - 1 



8N 



N - 1 



(31) 



Figure 9 shows a plot of the relative error of the predicted values of 
D versus log 8 2 for fix = 0.01, 0.005, 0.002, • • • , 0.00002. Except for 
a small region on the curve <5i = 0.01, the relative error is less than 1.0 
percent for the entire range of <5 2 and <5i considered. Based on a value of 
N = 127, a relative error of 1.0 percent in D is equivalent to an error 
of 1.26 samples in N, or approximately one-half a filter order off from 
the correct order. Errors of this magnitude are generally considered to 
be quite small, i.e., the prediction is reasonably good. 

In an effort to improve the fit and extend the range of applicability 
of the approximation, a nonlinear formula was chosen for D. Based on 
the data of Fig. 8, it was observed that the slope of the curve of D 
versus logio 5 2 changes nonlinearly with log fii. The simplest approxi- 
mation was to try a fit which was linear with respect to log 8 2 and 
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quadratic in log Si. Such a fit is of the form 

#nl = [ai(logio <5i) 2 + «2 login 5i + 03] logio 6 2 

+ [o 4 (logio «i) 2 + a 5 log M Si + a fl ]. (32) 

The constants ai to a& were chosen to minimize the mean-square rela- 
tive error for the N = 127 data over the range 0.1 ^ Si ^ 0.000001, 
0.1 ^ <5 2 ^ 0.000001, and turned out to be 

ai = 5.309 X 10~ 3 
a 2 = 7.114 X 10" 2 
a 3 = -4.761 X 10" 1 
a 4 = -2.660 X 10" 3 
a 5 = -5.941 X 10" 1 
a 6 = -4.278 X 10" 1 . 

Figure 10 shows the relative error of the predicted value of D as a func- 
tion of log 8 2 for values of 5i from 0.1 to 0.00002. The peak percentage 
error is 1.3 percent, and over most of the range the percentage error is 
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Fig. 9 — The relative errors in fitting the data of Fig. 8 with a linear curve over the 
range 5i ^ 0.01, for various values of 5i (N = 127). 
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Fig. 10 — The relative errors in fitting the data of Fig. 8 with a nonlinear curve over 
the range 5i ^ 0.1, for various values of 5i (N = 127). 



much smaller. Clearly, this prediction formula is acceptable for almost 
any design application. 

Figure 11 shows a summary of the predicted values of D as a func- 
tion of stopband attenuation for a wide range of values of passband 
ripple. In this case, passband ripple in dB is defined as 



where 



and 



Passband ripple = 20 logio f , _ - j , 
Stopband attenuation = — 20 logio (8.), 
2«i 



5 B = 



5 a = 



1 + Si 

1 + *l" 



(33) 
(34) 

(35) 
(36) 



These data correspond to standard design data for continuous-time 
filters where the frequency response magnitude is constrained to be less 
than or equal to 1.0. 
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Fig. 11 — Plots of D versus 201ogi 5 2 for various values of the parameter 
20 logio (1/1 — S p ), where S p = 25i/(l + 50, as calculated from the optimum non- 
linear fits to the data of Fig. 8 (N = 127). 



V. CORRECTIONS FOR SMALL VALUES OF N 

The formulas in the previous section are accurate for predicting D 
(or, equivalently, N) for values of N greater than about 51. As seen from 
the curves in Fig. 7, as JV decreases, D decreases for fixed values of 52 
and K. It is also seen from Fig. 7 that the differences increase with 
decreasing <5 2 or, equivalently, increasing transition width. An examina- 
tion of the relative deviation of Z> NL , the predicted value of D, from its 
true value as a function of transition bandwidth showed that, inde- 
pendent of N, the deviations could be simply approximated by a curve 
of the form 

D COT = f(K)(AF)>, (37) 



where D C0T is the correction term and f(K) is of the form 
f(K) = (0.51244 log 10 K + 11.01217). 



(38) 



(The constants in eq. (38) were again obtained by a minimum mean- 
square relative error data-fitting procedure.) Thus, using eqs. (37) and 
(38), a formula for D which depends on N, K, and 8 2 can be obtained. 
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Adopting the notation 

(i) D m (8 2 ,K) = D NL of eq. (32)* 
(ii) t>{8 2 , K, N) = predicted value of D as a function of N as well 

as 8 2 and K 
(in) D(8a, K, N) = true value of D, 

we obtain the relation 

D„(h, K) - D{8 2 , K, N) = f(K)(AF)\ (39) 

Thus, in a design case where h, K (or 81), F p , and F, (i.e., AF) are 
specified and the problem is to estimate N, the impulse response dura- 
tion required to meet these specifications, eq. (39) can be used directly 

since 

D(8 2 ,K,N) = (N - 1)AF. (40) 

Thus, combining eqs. (37) and (38) and solving for N gives 

N _ DM2, K) - f(K)(AFy , . 

iV AF 

A closed-form expression for £>(8 2 , K, N) may be obtained by sub- 
stituting eq. (40) for £>(8 2 , K, N) in eq. (39) and solving the quadratic 
equation for AF. Equation (40) is then used to give 2)(5 2 , K, N). Thus, 

we get , v 

f(K)(AF) 2 + (N - 1)AF - DM2, K) = (42) 

In the limit, as N tends to infinity, eq. (44) shows t){8 2 , K, N) tends to 
D K (8 2 , K) as expected. 

Using eq. (44), the relative error of the predicted value of D from the 
true value was measured for the data for N = 3, 7, 11, 19, 51, and 127 
with 81 ^ 0.1. The relative errors for values of K from 1 to 500 are 
plotted in Figs. 12a through 12f. In all cases, the worst relative error in 
D is sufficiently small that the equivalent error in N is less than one 
sample. Thus, for all practical purposes, the design equations above 
serve as a useful guide for estimating the order of the filter required to 
meet design specifications. 

* Eq. (32) gives Dnl as a function of 5i and 5 2 but since &i = K5 2 , it is also im- 
plicitly a function of 5 2 and K. 
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VI. SUMMARY OF DESIGN PROCEDURES AND EXAMPLES 

Based on the results presented in this paper, it is now possible to give 
a set of rules for estimating the filter impulse response duration re- 
quired to meet given design specifications. These rules are as follows : 

(i) Check if either «5i or 82 is greater than 0.1, in which case the 

graphical data of Fig. 7 are used directly to estimate N. 
(ii) Calculate a value of N, call this Ni, corresponding to the extra- 
ripple case where N p = N e , from the equation 

Ni ~ AF + l ' 

where D M (di, 8 2 ) is the optimum nonlinear fit to the data for 
large N and is given by : 

Ao(*i,«a) = [5.309 X 10" 3 [logio 8iJ + 7.114 X lO-Mogio^ 
- 0.4761] log 10 5 2 - [2.66 X lO" 3 (log 10 5i) 2 

4- 0.5941 log 10 S l 4- 0.4278] 
and f(K) is given by 

f(K) = 0.51244 log 10 (K) + 11.01217. 

(Hi) If the desired value of F p is less than or equal to 0.04 (let us 
call this case 1), or if the desired value of F a is greater than or 
equal to 0.46 (case 2), then the estimate of N is obtained in 
rule (iv). Otherwise, the value Ni of rule (ii) is used as the 
estimate of N. 

(iv) To obtain the value of N c for the Chebyshev solution for case 
1, eq. (26) is used to get a first approximation to N c . N c is then 
systematically varied until a Chebyshev solution is obtained 
which meets specifications on $1, 8 2 , and AF. (As discussed 
earlier, it is not generally possible to find a Chebj'shev solution 
which meets specifications exactly on all four filter parameters.) 
For case 2, 81 and 82 are interchanged, and F a is replaced by 
0.5 — F p , in order to solve for N c from eq. (26). In cases where 
this rule is applied, the value N c obtained is a lower bound to 
the true value of Ni. No upper bound may be given in this case. 
A discussion of this problem is given below. 

Several comments are necessary about these rules before proceeding 
to some examples. The discussion in this paper has concentrated on two 
regions of the curve of AF versus F p - the Chebyshev solutions and the 
case of extraripple filters with an equal number of passband and stop- 
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band ripples. The j ustification for such emphasis was the results shown 
in Figs. 2 and 3 which indicated that the extraripple solution for 
N p = N. was fairly representative of a large region of the curve of 
AF versus F p , except in the case where F p was either very small or very 
large, in which case the Chebyshev solution became important. Also, 
as seen from Fig. 2, between extraripple solutions, the curve of AF 
versus F p peaks up. However, it is seen that in many cases one can 
"approximately" bound the maximum between extraripple solutions 
by the next lower-degree extraripple solution. Since the design equa- 
tions can predict AF for this case (the next lower-degree extraripple 
solution) , a good bound on N can be obtained for a reasonably large 
region of the curve of AF versus F p . In these cases, the value given by 
rule (it) is good to within ±4 in the worst cases, i.e., large values of 
K, and generally to within ±2. 




PASSBAND CUTOFF FREQUENCY IF 



Fig. 13— Explanation of the types of behavior of curves of transition width versus 
passband cutoff frequency. 



OPTIMUM FIR LOW-PASS DIGITAL FILTERS 



795 



In the region of the Chebyshev solutions, however, such bounding 
procedures no longer are valid. An explanation of the difficulties which 
may be encountered is given in Fig. 13, which shows two types of curves 
of AF versus F p . Fig. 13a shows the case, discussed above, where the 
extraripple solutions of impulse response duration N — 2 are approxi- 
mately midway between extraripple solutions of impulse response dura- 
tion N. In this case, the next degree solution bounds the maximum AF 
between extraripple solutions. Figure 13b shows the case where the 
extraripple solutions of impulse response duration JV — 2 have ap- 
proximately the same values for AF and F p as the extraripple solutions 
of impulse response duration N. In these cases, there is no good bound 
on the maximum value of AF between adjacent extraripple solutions. 
The case of Fig. 13b corresponds to regions of F p near 0.0 and F, near 
0.5, i.e., in the regions of the Chebyshev solutions. In these cases, as 
discussed in rule (iv), there is only an underbound on N, and no 
overbound. 

The choice of a value of 0.04 in rule (in) as the width of the region 
during which the behavior of N can only be underbounded was ob- 
tained from the data of Fig. 3 which shows that, beyond this region, the 
variation in the values of AF for extraripple solutions is small. 

Figures 14 through 16 illustrate typical behavior of the curve of 
minimum value of impulse response duration N, to meet given specifi- 
cation on Si, 82, and AF as a function of F p . Figure 14 shows data for 
the case 5 X = 0.01, 8 2 = 0.0001, AF = 0.158. The value of Ni from 
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Fig. 14 — Optimum values of N to meet given design specifications on 5i, 5 2 , and 
AF, as a function of F p . 
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rule (ii) is 21, and the values of N e are 19 (for small F p ) and 13 (for 
large F p ). Over the region 0.06 ^ F p g 0.32, N is within 2 of the nomi- 
nal value of N\ = 21. Over the entire range of F p , N is within 4 of the 
value 2Vi = 21. The data of Fig. 14 correspond to the case of Fig. 13a. 
Figure 15 shows data for the case 5 X = 0.01, S 2 = 0.0001, AF = 0.032. 
The value of Ni from rule (ii) is 101, and the values of N e are 95 (for 
small F p ) and 55 (for large F p ). Over the region 0.1 g F p g 0.38, N is 
within 2 of the nominal value of iV\ = 101. However, in the region 
^ F p ^ 0.036, the value of N fluctuates from a minimum of 95 (the 
Chebyshev lower bound) to a maximum of 117. The explanation of this 
erratic behavior of N is seen in Fig. 16. Figure 16a shows the data of 
Fig. 15 on an expanded horizontal scale, and Fig. 16b shows a plot of 
the approximate curves of AF versus F p for all values of N from 95 to 
117. The heavily traced parts of these curves show the lowest-order 
solution which just meets specifications on AF. From Fig. 16b, it is 
clear that in the vicinity of the first few extraripple solutions, the curves 
of A^ versus F p are exceedingly steep. Hence, a slight change in F p 
greatly increases the required order solution to meet specifications. In 
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Fig. 16— Explanation of the behavior of the data of Fig. 15 in the region ^ F r 
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these cases, it is impossible to estimate the exact filter impulse response 
duration which is required. Instead, only a lower bound can be given. 
Fortunately, as seen in Fig. 15, the regions in which this erratic be- 
havior can occur are limited. 

We conclude this section with a set of examples which illustrate the 
use of the design rules. 
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Example 1: Find the minimum value of N required to meet the specifi- 
cations 5! = 0.05, 8 2 = 0.0001, F p = 0.19, F, = 0.21. 

From rule (ii) we get Ni = 129.7 which is rounded to 129. (Herein, 
all values of N will be rounded to the nearest odd integer.) Since F p 
is far from the range for the Chebyshev solutions, the value of 129 is 
used as the appropriate estimate. The actual filter impulse response 
duration required is N = 131, although the N = 129 filter just missed 
meeting specifications. 

Example 2: Find the minimum value of N required to meet the specifi- 
cations Sx = 0.01, 8 2 = 0.0001, F p = 0.213, F, = 0.373. 

From rule (ii) we get 2Vi = 19. Since F p is again out of the range of 
the Chebyshev solutions, the value 19 is used as the estimate of N. 
The actual filter impulse response duration required is AT = 19. 

Example 3: Find the minimum value of N required to meet the specifi- 
cations 5i = 0.1, 8 2 = 0.1, F p = 0.12, F, = 0.19. 

From rule (ii) we get N 1 = 11. Since F p is out of the range of the 
Chebyshev solutions, the value 11 is used as the estimate of N. The 
actual value of AT is 11. In this case, it is interesting to note that the 
value of N e for the Chebyshev solution is also 11. This example points 
out that, for filter specifications leading to small values of N, there is 
very little variation in the actual value ofN&sF p varies. This observa- 
tion has been made earlier with respect to Fig. 14. 

Example 4: Find the minimum value of N required to meet the speci- 
fications 5i = 0.01, 8 2 = 0.0001, F p = 0.36, F. = 0.497. 

From rule (ii) we get Ni = 23. Since F, is within the bounds of rule 
(Hi) (case 2), N c is computed from rule (iv) as 13. The actual value of 
AT is 19. In this case, the lower bound is within three filter orders of the 
true solution. 

VII. SUMMARY 

This paper has presented a wide variety of data on the relationships 
between design parameters for optimal low-pass FIR linear-phase digi- 
tal filters. Analytical formulas were derived for the Chebyshev solu- 
tions, i.e., when there was only one passband or stopband ripple. Ap- 
proximate fits to the data using nonlinear relationships between &F, 
N, 8 h and 8 2 were given in the case where the number of passband and 
stopband ripples was equal, and it was argued that these relationships 
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were valid over a wide range of values of the filter parameters. Finally, 
a simple set of rules for estimating the minimum value of N, which 
meets given specifications on F p , F 8 , 81, and 8 2 , was discussed. Examples 
were given to illustrate the application of the rules. 
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